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Spin systems with long-range interactions are "non-extensive" if the strength of the interactions 
falls off sufficiently slowly with distance. It has been conjectured for ferromagnets, and more recently 
for spin glasses, that, everywhere in the non-extensive regime, the free energy is exactly equal to that 
for the infinite range model in which the characteristic strength of the interaction is independent of 
distance. 

In this paper we present the results of Monte Carlo simulations of the one-dimensional long-range 
spin glasses in the non-extensive regime. Using finite-size scaling, our results for the transition 
temperatures are consistent with this prediction. We also propose, and provide numerical evidence 
for, an analogous result for diluted long-range spin glasses in which the coordination number is finite, 
namely that the transition temperature throughout the non-extensive regime is equal to that of the 
infinite-range model known as the Viana-Bray model. 

PACS numbers: 05.50.+q 75.50.Lk 75.40.Mg 



I. INTRODUCTION 

In the theory of phase transitions, it is often helpful 
to study models in a range of dimensions from above the 
"upper critical dimension", d u , where mean- field critical 
behavior is valid, to below the "lower critical dimension" , 
d/, where fluctuations destroy the transition. For Ising 
spin glasses d\ ~ 2.5 |l| and d u = 6. However, it has 
been difficult to cover this broad range numerically for 
spin glasses, since d u is quite large, and slow dynam- 
ics prevents equilibration at low temperatures when the 
number of spins TV (= L d ) is greater than a few thousand. 
It follows that at and above d n , one cannot study a suf- 
ficient range of linear sizes L to perform the necessary 
finite-size scaling (FSS) analysis. 

As a result, there has been a lot of recent attention on 
long-range models in one dimension, in which the inter- 
actions fall off with a power of the distance. Such models 
have a long history going back to Dyson 0, Q , who con- 
sidered a ferromagnet with interactions Jij falling off like 
l/r^ 7 , and found a paramagnet-ferromagnet transition for 
a < 2. Kotliar et al. [4] were the first to study the spin 
glass version of this model, which has received a lot of 
attention numerically in the last few years [E|-[9|. 

Varying the power a in the long-range spin glass model, 
one has a range of behavior similar that obtained by vary- 
ing the dimension in short-range models, namely there is 
a "lower critical value", o\ = 1, [To| above which there is 
no transition at finite temperature, and an "upper criti- 
cal value", <j u = 2/3 [4], below which the transition has 
mean field critical exponents. Note that increasing a 
makes the interactions more short range, and so corre- 
sponds to decreasing d. 

A precise connection between d for short-range models 
and c for long-range models can be made in the mean 
field region (d > 6 or 1/2 < a < 2/3), namely [ll| 

(d>6, 1/2 <cr< 2/3). (1) 



This mapping shows that d— ^oofora— ^1/2. Since the 
transition temperature in mean field theory is given by 
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we see that for smaller values of a, i.e. < a < 1/2, 
the strength of the interactions has to be scaled with an 
inverse power of the system size to obtain a sensible ther- 
modynamic limit. We call this regime "non-extensive". 
The extreme limit of this region, a = 0, is the Sherrington 
Kirkpatrick (SK) model [12j, which is "infinite-range". 
To complete the picture of the 1-d long-range spin glass 
model, in this paper we study the non-extensive regime 
(0 < a < 1/2), which has not been studied before, to our 
knowledge, apart from the SK model (a = 0). 

The non-extensive regime for ferromagnets has already 
been investigated [l3|, Q • This work shows that the be- 
havior in the whole non-extensive regime is the same, 
with a suitable rescaling of the interactions, as that of 
the infinite-range ferromagnet in which every spin inter- 
acts equally with every other spin, i.e. a = 0. We give 
intuitive arguments for this in Appendix \K\ 

It is interesting to ask if the same is true for spin 
glasses. In a recent paper Mori [15] has claimed that 
this is so, i.e. for all < a < 1/2 the behavior is the 
same as that of the SK model (a = 0) provided the inter- 
actions are scaled with system size so that S^Jd^-Jav 
is set to the same value for all a. However, this argu- 
ment is just at the level of replicating the Hamiltonian 
so it becomes translationally invariant, and then arguing 
that the earlier work for ferromagnets can be taken over 
directly to prove the result. While plausible, this result 
is by no means rigorous and so we test it here by Monte 
Carlo simulations. 

One of the models we simulate here is the usual one in 
which every spin interacts with every other spin. How- 
ever, it is also interesting to carry out the same study 
for a diluted model [3] with a fixed average coordination 
number z. This model has received a lot of attention re- 
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cently because the computer time per sweep only varies 
as N z (rather than N 2 for the undiluted model), so it 
can be simulated much more efficiently than the undi- 
luted model for large N. The diluted model with a = 
is called the Viana-Bray [16| model. It corresponds to a 
spin glass on a random graph, the exact solution of which 
is expected to be the Bethe-Peierls approximation. By 
analogy with Mori's proposal, we suggest here that the 
behavior of the diluted spin glass model is identical to 
that of the Viana-Bray model a = everywhere in the 
non-extensive region (0 < a < 1/2). We shall also pro- 
vide numerical evidence for this. 

We should emphasize that universal quantities, such 
as critical exponents, are expected to be the same ev- 
erywhere both in the mean field (1/2 < a < 2/3) and 
non-extensive (0 < a < 1/2) regimes. The claim that 
we test is that all the behavior of these models (not just 
the critical behavior) is identical for all a in the non- 
extensive regime, at least in the thermodynamic limit. 
We therefore need to look at non-universal quantities, 
and focus here on one particularly convenient quantity, 
the value of the transition temperature T c . 

The plan of this paper is a follows: In Sec. |TT] we 
describe the models used in the simulations and give 
their corresponding mean-field transition temperatures. 
In Sec. IIIII we give the details of the Monte Carlo simula- 
tions and FSS analysis. The results are given in Sec. HVl 
and our conclusions are summarized in Sec.|Vl Appendix 
\K\ provides an intuitive explanation of why the behav- 
ior of the ferromagnet is independent of a in the non- 
extensive regime. 



II. MODELS 



The Hamiltonian that we study is 



H 



J ij $i Sj , 



(3) 



where the 5^, (i = 1,2,- •• , L) are Ising spins which 
take values ±1, and the are statistically independent, 
quenched, random variables. The mean is taken to be 
zero and the variance satisfies 
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where, for the distance we put the sites on a ring and 
take the chord distance between sites % and j [17j, i.e. 
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The form of the distribution of the Jij is different for 
the undiluted and diluted models. For the undiluted case 
the distribution of the is Gaussian, 



where the variance is given by 
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in which C is a constant to be determined below. 

In order to compare models with different values of <r, 
for each a and L, we scale the variance so that 
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(undiluted) , 
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where the sum is for fixed i and we have Ju = 0. Equa- 
tion ([8]) determines the value of C in Eq. ([7]). Because 

we consider the non-extensive regime, C must vanish for 

_ i 

L — >> oo like L a 2 . 

The expression for the mean-field transition temper- 
ature in Eq. (j2j) is the exact result for the SK model, 
cr = 0. Hence, from Eq. J5J), we have 



T c (cr = 0) = 1 (undiluted) 



(9) 



For the diluted model, rather than the strength of the 
interaction falling off like 1/rfj, most bonds are absent 
and it is the probability of there being a non-zero bond 
which falls of with distance (asymptotically like 1/rfj 7 ). If 
a bond is present it is chosen from a Gaussian distribution 
with mean zero and variance unity (i.e. independent of 
Vij). In other words 
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(diluted), 
(10) 

where pij ex l/rfj 7 at large distance. 

It is convenient to fix the mean number of neighbors z. 
The pairs of sites with non-zero bonds are then generated 
as follows. Pick a site i at random. Then pick a site j 
with probability = A/rfj 7 , where A is determined by 
normalization. If there is already a bond between i and j 
repeat until a pair j is selected which does not already 
have a bond 1 . At that point set equal to a Gaussian 
random variable with zero mean and variance unity. This 
process is repeated Nz/2 times so the number of sites 
connected to a given site has a Poisson distribution with 
mean z. Because each site has, on average, z neighbors, 
and the variance of each interaction is unity, we have 
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The transition temperature for the diluted model with 
cr = was shown by Viana and Bray [l6| to be given by 
the solution of 
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, (undiluted), 1 Note that if zpij < 1 then in Eq. JTU) is given by = zpij, 

but otherwise there are corrections due to rejection of pairs i,j 
(6) when there is already a bond between them. 
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L 


^Vsamp 


iVequil 


^Vmeas 




Trnax 


N T 





64 


16000 


1000 


10000 


0.5 


1.65 


47 





128 


16000 


1000 


10000 


0.5 


1.6 


45 





256 


16000 


1000 


10000 


0.5 


1.6 


45 





512 


8000 


1000 


10000 


0.75 


1.55 


33 





1024 


8000 


1000 


10000 


0.75 


1.5 


31 





2048 


4000 


1000 


10000 


0.75 


1.5 


31 





4096 


4000 


2000 


10000 


0.85 


1.525 


28 


0.25 


64 


16000 


1000 


10000 


0.5 


1.65 


47 


0.25 


128 


16000 


1000 


10000 


0.5 


1.6 


45 


0.25 


256 


16000 


1000 


10000 


0.5 


1.6 


45 


0.25 


512 


8000 


1000 


10000 


0.5 


1.525 


42 


0.25 


1024 


8000 


1000 


10000 


0.75 


1.5 


31 


0.25 


2048 


4000 


1000 


10000 


0.75 


1.5 


31 


0.25 


4096 


4000 


2000 


10000 


0.85 


1.525 


28 



TABLE I: Simulation parameters for the undiluted models. 
iVsamp is the number of samples, iVequii and N m eas are the 
number of sweeps for equilibration and for the measurement 
phase, respectively. We simulate Nt temperatures between 
^min and T max . 
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L 


-/Vsamp 


^Vequil 
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Tmin 


Tmax 


N T 





256 


8000 


400 


8000 


1.85 


2.5 


27 





512 


8000 


800 


16000 


1.85 


2.5 


27 





1024 


8000 


2000 


40000 


1.85 


2.5 


27 





2048 


4000 


2000 


40000 


1.85 


2.5 


27 





4096 


4000 


2000 


40000 


1.9 


2.5 


25 





8192 


2000 


4000 


80000 


1.9 


2.5 


25 





16384 


2000 


4000 


80000 


2.0 


2.5 


14 


0.25 


256 


8000 


800 


16000 


1.85 


2.5 


27 


0.25 


512 


8000 


800 


16000 


1.85 


2.5 


27 


0.25 


1024 


8000 


1200 


24000 


1.85 


2.5 


27 


0.25 


2048 


4000 


2000 


40000 


1.85 


2.5 


27 


0.25 


4096 


4000 


2000 


40000 


1.9 


2.5 


25 


0.25 


8192 


2000 


4000 


80000 


1.9 


2.5 


25 


0.25 


16384 


2000 


4000 


80000 


2.0 


2.5 


14 


0.375 


256 


32000 


1200 


24000 


1.863 


4.0 


24 


0.375 


512 


26327 


1200 


24000 


1.863 


4.0 


26 


0.375 


1024 


16000 


1200 


24000 


1.913 


4.0 


24 


0.375 


2048 


15998 


2000 


40000 


1.95 


4.0 


24 


0.375 


4096 


8000 


4000 


80000 


1.962 


4.0 


28 


0.375 


8192 


7999 


4000 


80000 


1.975 


4.0 


34 


0.375 


16384 


4000 


4000 


80000 


2.0 


2.51 


18 



TABLE II: Simulation parameters for the diluted models, 
The parameters are the same as in Table [J 



We choose z — 6 for which we find 

T c (cr = 0) = 2.0564 , (diluted) . (13) 



III. METHOD 

We perform Monte Carlo simulations on the models 
described in Sec. [Ill To speed up equilibration we use 
the parallel tempering (exchange) Monte Carlo method 



[18|. In this approach one simulates Nt copies of the 
spins with the same interactions, each at a different tem- 
perature between a minimum value T m i n and a maximum 
value T max . In addition to the usual single spin- flip moves 
for each copy, we perform global moves in which we in- 
terchange the temperatures of two copies at neighboring 
temperatures with a probability which satisfies the de- 
tailed balance condition. In this way, the temperature 
of a particular copy performs a random walk between 
T m i n and T max , thus helping to overcome the free energy 
barriers found in the simulation of glassy systems. 

For the simulations of the undiluted model to be in 
equilibrium the following equality must be satisfied [l^ , 

(rpMF\ 2 

U= 2T (1 ~ m) ' ( undiluted )> ( 14 ) 
where 

^ = -E^^)]av 5 (15) 

(iJ) 

is the average energy, and q\ is the "link overlap" defined 

by 




(undiluted). (16) 



in which T C MF is given by Eq. (j2]) (and here set equal 
to unity by the scaling of the interactions, see Eq. J5J)). 
Equation (fH|) is obtained by integrating by parts with re- 
spect to the Jij the expression for the average energy, and 
noting that the distribution is Gaussian. This equation 
is useful because, very plausibly, the two sides approach 
their common value from opposite directions [17j, so, if 
the two sides agree, the system has reached equilibrium 
(at least for the energy and link overlap). 

For the diluted model, the equilibration test takes the 
form, [E[ 

U=-^(l-qi), (diluted), (17) 
where the link overlap is now defined by 

* = Wz ^ M^-> 2 )] av , (diluted), (18) 

(ho) 

in which = 1 if there is a bond between i and j, and 
zero otherwise. As with Eq. (fT4|) . we expect that the 
two sides of Eq. ([T7|) approach each other from opposite 
directions as equilibrium is approached. 

We consider results obtained by successively doubling 
the number of sweeps, in each case averaging over the 
last half of the sweeps, and we accept the data as being 
in equilibrium if the last three data points agree with 
each other within the error bars. The total number of 
sweeps used in this check is shown as iV equ ii in Tables [I] 
and [III We then do "production" runs where, in addition 



to A^equii sweeps for equilibration, we do 10 to 20 times 
as many sweeps, iV meas , during which measurements are 
performed. All the parameters used in the simulations 
are given in Tables HI and ITTI To avoid bias, each distinct 
thermal average, for example in Eq. ([T6]) . is evaluated 
in a separate copy (replica) of the system with the same 
interactions. 

We focus on moments of the spin glass order parameter 
q where 



t5> 



(1) o(2) 



(19) 



in which "(1)" and "(2)" refer to two independent copies 
of the system with the same interactions. Of particular 
interest are the spin glass susceptibility 



Xsg = L(q 2 ) , 



and the Binder ratio 
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Since the Binder ratio is dimensionless, its finite-size 
scaling (FSS) behavior is simple. We are always in the 
regime of mean field critical exponents (0 < a < 2/3), so 
it has the form [l9| 



g = g((T-T c )L 1 



/3 



(22) 



The spin glass susceptibility is not dimensionless but, 
since we are in the mean field regime, its FSS form is 
also known exactly. It has the form [19| 



X sg = L 1 ' z x({T-T c )L 1 ' z ) 



(23) 



One can therefore determine the transition tempera- 
ture from where the data for g or xsg/L 1 ^ 3 for different 
sizes intersects. However, we shall see that the data does 
not all intersect at a single temperature, showing that 
there are corrections to the FSS form in Eqs. ([22]) and 
([23]) . Consider Eq. ([23]) . According to standard finite- 
size scaling, the spin-glass susceptibility normally varies 
near the critical point according to 

XSG& L) = L a [f(L b t) + L~ u g(L y t) + • • • ] +co+d*+. • • : 

(24) 

where t = T — T c , a = 2 — n (= 2a — 1 here), and b = 
\jv. The L~ u term is the leading singular correction to 
scaling and Co is the leading analytic correction to scaling. 
However, in the mean field limit, a < 2/3, the exponents 
a and b are independent of aj20l-[22| and take the value at 
a c for all 1/2 < a < cr c , i.e. a = b = 1/3. Furthermore, 
although the L 2cr_1 term is replaced as the largest term 
by an L 1 / 3 term (due to the presence of a "dangerous 
irrelevant variable," cf. Refs. |20l-[22|) we expect [11[ this 



term to not disappear but rather become a correction to 
scaling. Hence, we replace Eq. ([M]) by 

Xsafa L) = L 1 / 3 [f{L^H) + L- u g(L^H) + • • • " 

+ doL^^hgiL^H) + c + c\t • • • . (25) 

The correction exponent uo can be obtained in the mean- 
field regime from the work of Kotliar et al. [4j and is given 
by uo = 2 — 3a. Hence, in the non-extensive regime, a < 
1/2 the dominant correction to scaling is the constant cq. 

Adding a constant to the RHS of Eq. ([23]) it is straight- 
forward to show that the intersection temperature of the 
data for xsg/L 1 ^ 3 for sizes L and 2L is given by 



T*(L,2L) = T C 



A 



(26) 



where A is a constant and the omitted terms are higher 
order in 1/L. We expect that the intersection tempera- 
tures for the data for g have the same form. We shall use 
Eq. ([26]) to determine T c for the models studied. 



IV. RESULTS 

We first present our results for the undiluted model. 
Data for the the scaled spin glass susceptibility and the 
Binder ratio are shown in Fig. [U The top part is for the 
SK model, a = 0, and the bottom part is for the undi- 
luted model with a = 0.25. One sees large corrections 
to scaling for the Binder ratio (the left-hand figures) but 
much smaller corrections for the scaled spin glass suscep- 
tibility (the right-hand figures). The inset enlarges the 
region of the intersections for the latter data. 

Figure [2] shows values for the intersection tempera- 
ture. These were determined by interpolation using cubic 
splines, and error bars computed by a jackknife analysis. 
For both values of a the data extrapolates to a value of 
1, the exact value for the SK model, (with very small er- 
rors). The quality of the fit, as represented by the good- 
ness of fit parameter Q [23[, is satisfactory except for 
the Binder ratio data for the SK model. We don't have 
a good explanation for this, except perhaps that multi- 
ple corrections to scaling are significant for the range of 
sizes studied. In any case we note that the result T c = 1 
for the SK model is rigorously correct. The result that 
T c = 1 also for a = 0.25, at the midpoint of the non- 
extensive region, provides strong evidence for the claim 
of Mori [lH that all models in the non-extensive region 
are identical to the SK model. While it would be useful 
to check this also in the space glass phase below T c , such 
simulations would be difficult because relaxation times 
increase dramatically at low T and so the range of sizes 
that could be studied would be much more limited than 
in the data presented here. 

The corresponding results for the diluted model for 
cr = and 0.25 are shown in Figs. [3] and [H We also 
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FIG. 1: (Color online) Data for the the scaled spin glass susceptibility and the Binder ratio for the undiluted model with a = 0, 
the SK model, (top), and the undiluted model with a = 0.25 (bottom). The dashed vertical line shows the exact value of the 
transition temperature (T c = 1) for the SK model. 



performed simulations for a = 0.375 and show the re- 
sulting intersection temperatures in Fig. [5] For a = 0, 
the Viana-Bray [l6| model, the transition temperature 
is given by Eq. (fl2j) which, for z = 6 taken here, gives 
the result in Eq. (fT3j) . In Fig. [3] we again see that cor- 
rections to scaling are larger for the Binder ratio than 
for the scaled spin glass susceptibility. The intersection 
temperatures all extrapolate to the exact value for a = 



within statistical uncertainty 2 



2 All the results within one standard deviation except for the data 
for g for a = and XSg/L 1 / 3 for a = 0.25 but even these 



are within 
acceptable. 



1.5 standard deviations, which we also consider 



6 




V. SUMMARY AND CONCLUSIONS 

We have performed Monte Carlo simulations to investi- 
gate the transition temperature of one-dimensional Ising 
spin glasses, both undiluted and diluted, for several val- 
ues of a in the non-extensive regime < a < 1/2. For 
the undiluted model we studied two values of a, a = 
and a = 0.25. For a = 0.25, which lies in the middle 
of the non-extensive regime, we find that the transition 
temperature agrees to high precision with the exact so- 
lution of the SK model. As a check, we also simulated 
the a = case, obtaining results consistent with the ex- 
act SK model result, though there seem to be multiple 
corrections to FSS for some of the data. 

For the diluted model we studied three values of a: 
a = 0, which corresponds to the Viana-Bray model, 
a = 0.25, which lies in the middle of the non-extensive 
regime, and a = 0.375. In all cases we found the tran- 
sition temperature to be consistent with the exact solu- 
tion of the Viana-Bray model; all results are within ~ 1.5 
standard deviations. 

To conclude, our results provide confirmation of the 
proposal [15[ that the behavior of (undiluted) spin glasses 
everywhere in the non-extensive regime is identical to 
that of the SK model. We have also proposed that an 
analogous result applies to diluted spin glass models, and 
provided numerical evidence for this too. 
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Appendix A: Spherical Approximation for the 
ferromagnet 

For the infinite-range ferromagnet, the interaction Jij 
is equal to 1/(N — 1) for i ^ j and for i = j. This 
fixes T c = 1. The Fourier transform of this interaction is 
given by 

J(k) = 4,o - ^ , (Al) 

so only the k = mode contributes to the transition. 

For a power-law decay of the interactions in the non- 
extensive regime (0 < a < 1), on dimensional grounds 
there is a singular piece which diverges like k a ~ 1 for k — »• 
0. Furthermore the interactions have to be multiplied by 
a number of order TV* 7-1 in order to satisfy the condition 
jiMF _ ^ J { - = 1. Hence, roughly speaking, we have 

J(k) oc (kNf- 1 , (A2) 

where we note that k = k n = 27rn/L, (n = 0, 1, 2, • • • ). 
(For n = 0, J(0) does not actually diverge but will be 
comparable to J(k\)). Hence other long wavelength 
modes, in addition to k = 0, are now significant. How- 
ever, we shall now see that there are not enough of them 
to change the value of T c from that of T C MF (= 1). 

We will do this by considering the "spherical approx- 
imation" [24|. in which we reexpress the problem as a 
Gaussian one, with "soft" spins <pi which take values from 
— oo to oo, and a Hamiltonian 

^Gauss = \^2 ~ ^0 ' ( A3 ) 

where \i is a Lagrange multiplier whose value is chosen 
to enforce the length constraint 

<#> = !■ (A4) 
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FIG. 3: (Color online) Data for the the scaled spin glass susceptibility and the Binder ratio for the Viana-Bray model, i.e. the 
diluted model with a = (top), and the diluted model with a = 0.25 (bottom). The dashed vertical line shows the transition 
temperature for the Viana-Bray model obtained from Eq. (|12j) of the text. 



It turns out the the spherical approximation is exact for 
an m-component model in the limit m — >• oo [25[ . Fourier 
transforming Eq. (|A4j) and doing the Gaussian integrals 



gives 



- = -Y, 



L^v- J(k) 



(A5) 



The transition occurs when the denominator vanishes at 



k = 0, i.e. when \i = J(0) and so 



nspher 



J(0) - J(k) ' 



(A6) 



It is interesting to compare this with the mean field re- 
sult, T c = J2j Jij = J(0). Since J u = (1/L) ^ fe J(k) = 



8 




FIG. 4: (Color online) Results for the intersection temperatures for the Viana-Bray model (left) and the diluted model with 
a = 0.25 (right). 
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FIG. 5: (Color online) Results for the intersection tempera- 
tures for the diluted model with a = 0.375. 



with the MF result being exact for this model). 

We now estimate T c from the spherical approximation, 
Eq. (|A6|h for the power-law model, where J(k) varies 
like Eq. (jA2|) . Because we normalize the interactions to 
J(0) = 1, we can include an extra factor of J(0) and 
expand in powers of J(fc)/J(0), i.e. 



— = -T- 

pher £ i 
k 

4e 



i 



J(fc)/J(0) 
J(k) [J(k)Y 

j(o) + U(o); 



= i 



/J(fc)Y 

U(o)y 



3 l 3 



fc J .)' 



(A8a) 
(A8b) 
(A8c) 
(A8d) 



we can rewrite the mean field transition temperature as 

k 

Thus, whereas in mean field theory, T c is equal to the 
average of J(0) — J(fc), in the spherical approximation 
1/T C is equal to the average of the inverse of this. 

For the infinite range model, where J(k) is given by 
Eq. (|A1|) and only the /c = mode contributes, the spher- 
ical result agrees with the mean field result (consistent 



where in Eq. (|A8c|) we used that ^ k J(k) = 0. In 
Eq. (lA8d|) we have Ej^j «: L 1_cr while Ej^ 2 j = 
L 1_2cr (0 < a < 1/2), = const. (1/2 < a < 1). 

Hence the leading correction term in Eq. (|A8d|) vanishes 
everywhere in the non-extensive regime. 

To conclude, in this appendix we have given a sugges- 
tive argument as to why T c for the ferromagnet is given 
exactly by the mean field value everywhere in the non- 
extensive regime. It is therefore also plausible that other 
properties are also identical to those of mean field theory 
(i.e. the infinite-range model.) 
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